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SUMMARY 


The  interaction  between  turbulent  fluid  motion  and  finite-rate  chemical  reactions  has 
not  been  studied  under  conditions  typical  of  a  high-temperatme,  hypersonic  boundary 
layer.  These  flows  are  significantly  different  than  turbulent  combustion  flows,  which  have 
been  studied  extensively.  In  hypersonic  flows  the  dominant  chemical  reactions  are  the 
dissociation  and  recombination  of  nitrogen  and  oxygen  molecules.  These  reactions  are  dif¬ 
ferent  than  combustion  reactions  for  several  reasons.  First,  the  equilibrium  composition  of 
reacting  air  depends  strongly  on  the  temperature.  This  is  contrasted  with  a  combustion 
process,  where  the  equilibrium  composition  is  determined  by  the  initial  fuel-oxidizer  ratio. 
Another  important  featme  of  air  reactions  is  that  they  huve  a  large  activation  energy,  and 
therefore  the  reaction  rate  is  typically  temperature  limited.  In  this  situation,  the  reaction 
rate  depends  exponentially  on  the  temperature,  and  small  increases  in  the  temperature 
result  in  large  increases  in  the  reaction  rate.  This  is  contrasted  with  non-premixed  com¬ 
bustion  flows  where  the  fuel-oxidizer  mixing  rate  determines  the  rate  of  product  formation, 
and  the  reaction  process  is  relatively  insensitive  to  the  temperature.  In  premixed  combus¬ 
tion  flows,  the  flame  speed  is  primarily  determined  by  the  diffusive  transport  of  reactive 
radical  species  and  heat,  rather  than  by  the  reaction  rate  itself. 

The  surface  of  a  hypersonic  vehicle  is  typically  cooled  below  the  adiabatic  wall  tem¬ 
perature.  Thus,  the  maximum  temperature  in  a  hypersonic  boundary  layer  occurs  some 
distance  from  the  surface  where  there  is  significant  shear  heating  and  the  wall  cooling 
is  not  important.  In  this  region,  the  air  molecules  dissociate,  and  the  reaction  products 
diffuse  toward  the  surface  where  they  recombine  due  to  the  cool  wall.  Therefore,  hyper¬ 
sonic  boundary  layers,  unlike  combustion  flows,  have  regions  of  endothermic  (dissociation) 
reactions  and  exothermic  (recombination)  reactions. 

Thus,  the  purpose  of  this  project  was  study  how  turbulent  temperature  fluctuations 
interact  with  the  finite-rate  chemical  reactions  that  occur  in  hypersonic  boundary  layers. 
The  compressible  Navier-Stokes  equations,  that  have  been  extended  to  include  the  effects 
of  finite-rate  chemical  reactions  and  heat  release,  are  solved  numerically.  This  allows 
the  study  of  how  changes  in  the  reaction  rate,  heat  release,  and  turbulent  Mach  number 
affect  the  decay  of  isotropic,  homogeneous  turbulence.  Isotropic  turbulence  is  the  most 
fundamental  turbulent  flow,  and  it  is  a  useful  idealized  test  case.  It  also  serves  as  an 
approximation  of  the  small-scale  turbulent  motion  in  a  boundary  layer. 

This  report  summarizes  the  results  of  the  AFOSR-sponsored  project,  “Simulation  of 
Turbulent  Hypersonic  Flows.”  More  details  of  the  work  may  be  found  in  the  two  attached 
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papers.  The  first  paper,  “Effect  of  Chemical  Reactions  on  Decaying  Isotropic  Turbulence,” 
describes  the  results  of  a  series  of  direct  numerical  simulations  of  compressible,  turbulent, 
reacting  flows  that  was  performed  as  part  of  this  project.  The  second  paper,  “Subgrid-Scale 
Model  for  the  Temperature  Fluctuations  in  Reacting  Turbulence,”  discusses  the  use  of  the 
simulation  data  to  develop  a  model  for  turbulent  temperature  fluctuations  in  hypersonic 
boundary  layers.  This  model  is  useful  for  large-eddy  simulations  of  hypersonic,  turbulent 
boundary  layers. 

The  research  was  carried  out  by  the  Principal  Investigator,  Graham  V.  Candler,  and 
M.  Pino  Martin,  a  Graduate  Research  Assistant  in  the  Aerospace  Engineering  and  Me¬ 
chanics  Department  at  the  University  of  Minnesota.  Two  conference  publications  resulted 
from  this  work:  M.P.  Martin  and  G.V.  Candler,  “Effect  of  Chemical  Reactions  on  the 
Decay  of  Isotropic  Homogeneous  Turbulence,”  AIAA  Paper  No.  96-2060,  June  1996;  and 
M.P.  Martin  and  G.V.  Candler,  “Evaluation  of  a  Subgrid-Scale  Model  for  the  Temperature 
Fluctuations  in  Reacting  Turbulence,”  AIAA  Paper  No.  97-0751,  Jan.  1997.  Also,  this 
work  was  presented  at  the  1996  and  1997  American  Physical  Society  Division  of  Fluid 
Dynamics  Meetings.  These  two  papers  have  been  revised  and  submitted  to  Physics  of 
Fluids. 


KEY  FINDINGS 

The  direct  numerical  simulations  (DNS)  of  turbulent  reacting  flows  showed  that  in 
the  case  of  exothermic  (heat-producing)  reactions,  a  portion  of  the  chemical  energy  is 
placed  in  the  compressible  turbulent  kinetic  energy  modes.  Almost  none  of  this  energy 
is  placed  in  the  incompressible  turbulent  kinetic  energy  modes.  Therefore,  exothermic 
turbulent  flows  tend  to  have  a  very  large  compressible  energy  component,  which  results 
in  strong  compressions  and  expansions  in  the  flow.  When  the  reactions  are  endothermic 
(heat- absorbing) ,  energy  is  removed  from  the  compressible  energy  modes. 

The  reacting  flow  that  was  studied  was  a  simple  perfectly  mixed  three-dimensional 
periodic  box  of  gas  that  was  initialized  with  realistic  isotropic  turbulence.  The  turbulent 
flow  was  then  allowed  to  react  and  decay  over  time.  Because  this  is  a  perfectly  mixed 
flow,  there  is  no  turbulent  length  scale  that  is  introduced  into  the  reaction  process.  This 
causes  the  addition  and  removal  of  energy  by  exothermic  and  endothermic  reactions  to 
be  scale-independent.  This  accounts  for  the  energy  affecting  only  the  compressible  energy 
modes. 

The  strength  of  the  turbulence- chemistry  interaction  depends  on  the  amount  of  heat 
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released  (or  absorbed)  and  on  the  rate  of  the  reactions.  Generally,  the  faster  the  reaction 
rate,  the  stronger  the  interaction.  This  shows  that  the  chemistry-turbulence  interaction 
depends  on  localization  of  the  heat  release  and  a  self-strengthening  feedback  mechanism. 
For  example,  consider  a  region  in  the  flow  where  there  is  a  positive  temperature  fluctuation. 
If  the  reaction  is  relatively  fast,  the  heat  released  by  the  reaction  will  remain  in  the  vicinity 
of  the  original  hot  spot,  resulting  in  a  further  increase  in  the  temperature.  However  if  the 
reaction  is  slow,  the  fluid  will  have  moved  before  the  heat  is  released  and  the  feedback 
process  will  be  weakened  or  eliminated. 

This  discussion  also  clarifles  why  the  chemical  energy  affects  the  compressible  kinetic 
energy  modes.  The  temperature  fluctuations  generated  by  the  reaction  rate  fluctuations 
produce  localized  regions  of  high  and  low  pressure.  This  results  in  compressive  and  expan¬ 
sive  (compressible  or  dilatational)  motion  of  the  gas. 

The  results  are  discussed  in  much  greater  detail  in  the  attached  paper,  “Effect  of 
Chemical  Reactions  on  Decaying  Isotropic  Turbulence,”  which  has  been  submitted  to 
Physics  of  Fluids. 

The  results  of  the  direct  numerical  simulations  were  used  to  formulate  a  model  for 
turbulent  fluctuations  under  conditions  typical  of  a  hypersonic  cruise  vehicle.  This  subgrid- 
scale  model  is  designed  for  use  in  large-eddy  simulations  of  these  flows,  and  embodies  the 
key  features  of  the  chemistry-turbulence  feedback  mechanism  discussed  above. 

The  subgrid-scale  model  represents  the  temperature  fluctuations  with  a  Gaussian 
probability  density  function  (PDF)  whose  mean  and  variance  are  fitted  to  the  DNS  data 
over  the  conditions  of  interest.  This  is  valid  because  the  DNS  data  show  that  the  distribu¬ 
tion  of  T'/ (T)  (temperature  fluctuation  normalized  by  the  ensemble-averaged  temperature) 
is  Gaussian.  Also,  the  mean  of  this  distribution  is  zero,  which  means  that  only  the  variance 
of  the  distribution  must  be  fitted.  A  single  parameter,  X/is,  characterizes  the  strength  of 
the  interaction.  Here,  A  is  the  Taylor  microscale,  which  is  a  turbulent  length  scale;  and  £e 
is  the  acoustic,  expansion  length,  weighted  by  the  non-dimensional  reaction  rate.  For  the 
exothermic  reactions  tested  to  date,  the  temperature  fluctuation  PDF  can  be  characterized 
by  a  power-law  fit  of  the  form 


where  Ah°  is  the  non-dimensional  heat  relea^  (which  characterizes  the  strength  of  the 
reaction).  A  and  B  are  fitting  parameters,  which  depend  on  the  Arrhenius  reaction  rate 
parameters. 
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If  the  reaction  is  endothermic,  a  simpler  fit  may  be  used,  depending  solely  on  the 
turbulent  Mach  number,  independent  of  the  reaction  type. 

Preliminary  tests  of  this  model  for  several  air  reactions  under  a  wide  range  of  condi¬ 
tions  give  good  agreement  with  the  original  simulations.  Again,  the  details  of  the  model 
and  the  comparison  of  the  model  with  the  original  DNS  data  are  presented  in  the  attached 
paper,  “Subgrid-Scale  Model  for  the  Temperature  Fluctuations  in  Reacting  Turbulence.” 
Therefore,  this  model  has  the  potential  for  use  in  large-eddy  simulations  of  hypersonic 
boundary  layers. 
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There  are  many  studies  of  turbulent  combustion  flows,  however  there  is  no  analysis  of  the 
interaction  between  turbulent  motion  and  the  chemical  reactions  that  typically  occur  in 
hypersonic  flows.  In  this  case,  the  rate  of  product  formation  depends  almost  exclusively  on 
the  temperature,  and  small  temperature  fluctuations  may  produce  large  increases  in  prod¬ 
uct  formation.  To  study  this  process,  we  perform  direct  numerical  simulations  of  isotropic 
turbulence  decay  under  conditions  typical  of  a  hypersonic  flow.  We  flnd  a  positive  feedback 
mechanism  between  the  turbulence  and  the  exothermic  reactions.  In  high  temperature  re¬ 
gions  the  reaction  rate  increases,  further  increasing  the  heat  released  and  the  temperature. 
Simultaneously,  the  heat  released  increases  the  pressure,  causing  localized  expansions  and 
compressions  that  feed  the  turbulent  kinetic  energy.  The  Reynolds  stress  budget  shows 
that  the  feedback  occurs  through  the  pressure-strain  term.  The  strength  of  the  feedback 
depends  on  how  much  heat  is  released,  the  rate  at  which  it  is  released,  and  the  turbulent 
Mach  number.  The  feedback  is  negative  for  endothermic  reactions. 

I.  Introduction 

The  interaction  between  turbulent  fluid  motion  and  flnite-rate  chemical  reactions  has 
not  been  studied  under  conditions  typical  of  a  high-temperature,  hypersonic  boundary 
layer.  These  flows  are  signiflcantly  different  than  turbulent  combustion  flows,  which  have 
been  studied  extensively.  In  hypersonic  flows  the  dominant  chemical  reactions  are  the 
dissociation  and  recombination  of  nitrogen  and  oxygen  molecules.  These  reactions  are  dif¬ 
ferent  than  combustion  reactions  for  several  reasons.  First,  the  equilibrium  composition  of 
reacting  air  depends  strongly  on  the  temperature.  This  is  contrasted  with  a  combustion 
process,  vrhere  the  equilibrium  composition  is  determined  by  the  initial  fuel-oxidizer  ratio. 
Another  important  feature  of  air  reactions  is  that  they  have  a  large  activation  energ}',  and 
therefore  the  reaction  rate  is  typically  temperature  limited.  In  this  situation,  the  leaction 
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rate  depends  exponentially  on  the  temperature,  and  small  increases  in  the  temperature 
result  in  large  increases  in  the  reaction  rate.  This  is  contrasted  with  non-premixed  com¬ 
bustion  flows  where  the  fuel-oxidizer  mixing  rate  determines  the  rate  of  product  formation, 
and  the  reaction  process  is  relatively  insensitive  to  the  temperature.  In  premixed  combus¬ 
tion  flows,  the  flame  speed  is  primarily  determined  by  the  diff'usive  transport  of  reactive 
radical  species  and  heat,  rather  than  by  the  reaction  rate  itself.^ 

The  surface  of  a  hypersonic  vehicle  is  typically  cooled  below  the  adiabatic  wall  tem¬ 
perature.  Thus,  the  maximum  temperature  in  a  hypersonic  boundary  layer  occurs  some 
distance  from  the  surface  where  there  is  significant  shear  heating  and  the  wall  cooling 
is  not  important.  In  this  region,  the  air  molecules  dissociate,  and  the  reaction  products 
diffuse  toward  the  surface  where  they  recombine  due  to  the  cool  wall.  Therefore,  hyper¬ 
sonic  boundary  layers,  unlike  combustion  flows,  have  regions  of  endothermic  (dissociation) 
reactions  and  exothermic  (recombination)  reactions. 

Thus,  it  is  the  purpose  of  this  paper  to  study  how  turbulent  temperature  fluctuations 
interact  with  the  finite-rate  chemical  reactions  that  occur  in  hypersonic  boundary  layers. 
We  numerically  solve  the  compressible  Navier-Stokes  equations  that  have  been  extended  to 
include  the  effects  of  finite-rate  chemical  reactions  and  heat  release.  We  study  how  changes 
in  the  reaction  rate,  heat  release,  and  turbulent  Mach  number  affect  the  decay  of  isotropic, 
homogeneous  turbulence.  Isotropic  turbulence  is  the  most  fundamental  turbulent  flow, 
and  it  is  a  useful  idealized  test  case.  It  also  serves  as  an  approximation  of  the  small-scale 
turbulent  motion  in  a  boundary  layer. 

Many  direct  numerical  simulations  of  reacting  turbulent  flows  have  been  reported 
for  combustion  applications. However,  as  mentioned  above,  these  simulations  focus 
on  different  issues  than  the  current  paper.  For  example,  Gao  and  O’Brien®  and  Javery 
et  al}‘^  consider  the  interaction  of  chemical  reactions  with  homogeneous  turbulence,  but 
the  reaction  rate  is  taken  as  constant  and  the  effects  of  heat  release  are  not  considered. 
Picart  et  al}  consider  a  concentration-dependent  reaction  rate  to  simulate  the  exponential 
dependence  of  the  reaction  rate  on  temperature,  but  their  simulations  do  not  include  heat 
release  effects.  Several  simulations  of  reacting  shear  layers  have  been  performed  including 
the  effects  of  heat  release.  McMurtry  et  al}  and  Grinstein  and  Kailasanath®  use  a  constant 
reaction  rate,  Givi  et  al.‘  and  Miller  et  consider  both  constant  and  temperature- 

dependent  reaction  rates.  These  simulations  show  a  significant  effect  of  heat  release  on 
the  structure  of  the  mixing  layers  studied.  Therefore,  the  previous  simulations  have  some 
similarity  to  those  presented  here,  but  they  focus  on  different,  combustion-related  issues. 

In  the  remainder  of  the  paper,  we  introduce  the  equations  of  motion  and  relevant  non- 
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dimensional  parameters  for  a  mixture  of  reacting  gases.  We  then  discuss  the  formulation  of 
our  direct  numerical  simulations  and  present  the  results.  Finally,  we  analyze  the  chemistry- 
turbulence  interaction  process  and  summarize  our  findings. 


II.  Governing  equations 

The  equations  describing  the  unsteady  motion  of  a  reacting  flow  with  no  contribution 
of  vibrational  modes  are  given  by  the  species  mass,  mass-averaged  momentum,  and  total 
energy  conservation  ecjuations 

dps  ^  d  f  , 

dE  d  /  _ >  \ 

-^  +  ^  ((^  +  -  UiTij  +qj  +  Y^  PsV.jhsJ  =  0 

where  Ws  represents  the  rate  of  production  of  species  s  due  to  chemical  reactions;  pa  is 
the  density  of  species  s;  Uj  is  the  mass- averaged  velocity  in  the  j  direction;  Vsj  is  the 
diffusion  velocity  of  species  s;  p  is  the  pressure;  Tij  is  the  shear  stress  tensor  given  by  a 
linear  stress-strain  relationship;  qj  is  the  heat  flux  due  to  temperature  gradients;  ha  is  the 
specific  enthalpy  of  species  s;  and  E  is  the  total  energy  per  unit  volume  given  by 


£'  =  -t-  yuiUi  +  paK, 


(2) 


where  c^a  is  the  specific  heat  at  constant  volume  of  species  s;  and  h°  represents  the  heat 
of  formation  of  species  s. 

To  derive  the  expression  for  Wa,  consider  a  reaction  where  species  Si  reacts  to  form 
species  S2 

SI  +  M  ^  S2  -h  M,  (3) 

where  M  represents  a  collision  partner,  which  is  either  Si  or  S2  in  this  case.  The  source 
terms  for  Si  and  S2  can  be  written  using  the  law  of  mass  action 
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-)-  Msi  kb 
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and  WS2  =  — 'fusi-  kf  and  kf,  are  the  forward  and  backward  reaction  rates  respectively; 
These  are  written  as 


^  ’‘f 


(5) 


where  K^g  is  the  temperature-dependent  equilibrium  constant. 

For  a  two  species  mixture,  the  diffusion  velocity  can  be  accurately  represented  using 
Fick’s  law 


Ps'^sj  —  pD 


dcs 
dxj  ’ 


(6) 


where  =  p^fp  is  the  mass  fraction,  and  D  is  the  diffusion  coefficient  given  in  terms  of 
the  Lewis  number 

pDPr 


Le  = 


P 


(7) 


where  Pr  is  the  Prandtl  number,  p  is  the  viscosity,  and  Le  is  taken  to  be  unity,  so  that 
the  energy  transport  due  to  mass  diffusion  is  equal  to  the  energy  transport  due  to  thermal 
conduction. 


III.  Governing  parameters 

In  this  section,  we  introduce  the  non-dimensional  numbers  governing  the  physical  pro¬ 
cess.  Front  a  non-dimensional  analysis  of  the  governing  equations,  the  following  parameters 
are  obtained:  turbulent  Mach  number,  Mf,  Reynolds  number  based  on  the  lateral  Taylor 
microscale,  Rey,  Damkohler  number.  Da;  and  relative  heat  release,  Ah°.  The  turbulent 
state  of  the  flow  gives  rise  to  the  first  two  non-dimensional  numbers;  the  second  two  are 
present  if  the  flow  is  chemically  reacting.  These  parameters  may  be  written  as 


Da  = 


A|u;| 

pu'  ’ 


Rex 


pu'X 

? 

P 


Ah' 


CyT  -f- 
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where  q  is  the  RMS  magnitude  of  the  fluctuation  velocity;  a  is  the  speed  of  sound;  u'  is 
the  RMS  turbulent  velocity  fluctuation  in  one  direction;  Ah°  =  hg,  —  hgj  is  the  heat  of  the 
reaction;  and  c„  is  the  mixture  specific  heat  at  constant  volume. 

The  Damkohler  number  is  the  ratio  of  the  turbulent  time  scale,  ,  to  the  chemical  time 
scale,  Tc,  and  represents  a  non-dimensional  reaction  rate;  here  Tt  =  X/u'  and  Tc  =  pI\w\. 
Ah°  is  the  ratio  of  the  chemical  energy  released  or  absorbed  to  the  sum  of  the  internal  and 
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kinetic  energy.  A  negative  value  indicates  an  endothermic  reaction,  and  a  positive  value 
denotes  an  exothermic  reaction. 

IV.  Numerical  experiments 

It  is  the  purpose  of  this  paper  to  determine  how  turbulent  motion  interacts  with 
the  chemical  reactions  that  occur  in  hypersonic  boundary  layers.  A  logical  reaction  to 
consider  is  the  dissociation  and  recombination  of  nitrogen,  N2  +  M  ^  2N  +  M.  However, 
we  have  found  that  it  is  difficult  to  analyze  the  effects  of  this  reaction  on  the  decay  of 
isotropic  turbulence  because  many  parameters  change  simultaneously  during  the  reaction. 
For  example,  consider  a  case  where  we  start  with  pure  nitrogen  atoms,  and  let  them 
recombine  according  to  the  reaction  rate  and  equilibrium  constant  expressions  given  b}/ 
Gupta  et  Figure  1  plots  the  nitrogen  molecule  mass  fraction  for  simulations  with 

initial  temperatures  between  2000  and  7000  K  (corresponding  to  values  of  relative  heat 
release  between  33  and  9).  We  see  that  the  equilibrium  mass  fraction  depends  on  the 
initial  temperature  because  the  equilibrium  constant  is  a  function  of  temperature.  More 
importantly,  the  rate  of  reaction  differs  in  each  simulation.  The  differences  between  the 
simulations  make  it  difficult  to  compare  the  effects  of  turbulent  motion  on  the  evolution 
of  the  chemical  reactions.  Therefore,  we  must  devise  a  simplified  chemical  model  that 
captures  the  key  features  of  air  reactions,  yet  makes  it  possible  to  compare  simulations. 

There  are  several  issues  that  we  must  consider.  First,  if  the  average  molecular  weight 
of  the  gas  changes  when  the  reaction  occurs,  then  the  mixture- averaged  gas  constant  and 
specific  heats  change.  Therefore,  changes  in  the  pressure  and  internal  energy  will  not  only 
be  due  to  the  heat  released  by  the  reaction,  but  also  due  to  the  changing  properties  of 
the  gas  mixture.  To  circumvent  this  issue,  we  let  the  reactant  and  product  have  the  same 
molecular  weight  and  the  same  number  of  internal  degrees  of  freedom,  thus  the  mixture 
gas  constant  and  specific  heats  do  not  change  as  the  reaction  progresses.  For  the  same 
reason,  w^e  use  the  simple  binary  reaction  given  in  (3),  rather  than  a  dissociation  reaction 
in  which  the  average  molecular  weight  of  the  gas  changes  as  a  result  of  the  reaction. 

As  discussed  above,  the  equilibrium  constant  is  a  function  of  temperature,  which 
implies  a  different  chemical  equilibrium  state  for  simulations  with  diffeient  heat  release. 
Therefore,  we  take  Keg  to  be  constant,  so  that  the  equilibrium  chemical  state  is  fixed. 
This  allows  us  to  systematically  compare  simulations  with  different  initial  conditions  and 
values  of  heat  release. 

The  reaction  rate  is  also  a  function  of  temperature.  When  the  gas  temperature  goes  up 
the  reaction  rate  grows  exponentially,  which  presents  similar  problems  with  comparing  one 
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simulation  to  another.  Therefore,  we  must  modify  the  Arrhenius  form  of  the  reaction  rate 
given  by  (5).  It  is  desirable  that  the  modified  reaction  rate  be  constant  when  there  are  no 
temperature  fluctuations,  and  that  the  rate  vary  exponentially  when  there  are  temperature 
fluctuations.  One  such  reaction  rate  expression  is  given  by 

kfiT)  = 

6  T  If  e  T'\2  X  (^) 

~  ^  (T)  r  2  ( (T)  T  )  +•••)’ 

where  (T)  is  the  ensemble  average  of  the  temperature,  T'  =  T  —  (T)  is  the  temperature 
fluctuation,  C  is  a  constant,  and  9  is  the  activation  temperature  given  in  (5).  If  for  both 
reaction  rate  expressions,  we  compute  the  ratio  of  the  reaction  rate  at  temperature  T,  to 
that  at  the  ensemble- average  temperature,  (T),  we  obtain 


(10) 


Note  that  we  have  neglected  the  pre-exponential  factor  in  the  reaction  rate,  T'^,  because 
for  the  reactions  that  occur  in  hypersonic  flows,  ?;  <  —  1.  Thus,  at  large  temperatures,  this 
factor  has  a  weak,  sub-linear  variation  with  temperature.  Therefore,  the  two  rate  expres¬ 
sions  give  the  same  variation  of  the  relative  reaction  rate  for  any  temperature  fluctuation. 
The  only  difference  is  that  the  modified  reaction  rate  expression  gives  a  constant  rate 
when  the  temperature  is  equal  to  the  ensemble- average  temperature.  Thus,  the  modified 
reaction  rate  has  the  important  features  of  the  Arrhenius  rate  expression,  yet  allows  us  to 
systematically  compare  simulations. 


V.  Numerical  method  and  initial  conditions 


In  this  section  we  present  the  numerical  simulations  for  three-dimensional,  compress¬ 
ible,  homogeneous,  isotropic,  reacting  turbulent  flow.  The  numerical  method  used  is  based 
on  the  method  of  Lee  et  al}^  We  use  a  sixth-order  accurate  finite-difference  method 
based  on  a  compact  Fade  scheme  and  fourth-order  accurate  Runge-Kutta  time  integra¬ 
tion.  The  simulations  were  performed  on  grids  with  96^  points  following  the  criteria  given 
by  Reynolds^®  and  a  resolution  of  Krj  =  1  at  the  end  of  the  simulation,  where  K  is  the 
maximum  wavenumber  resolvable  on  the  grid  and  rj  is  the  Kolmogorov  scale.  The  compu¬ 
tational  domain  is  a  periodic  box  with  non-dimensional  length  'Itt  in  each  direction.  The 
velocity  field  is  initialized  to  an  isotropic  state  prescribed  by  the  following  energy  spectrum 


E(k)  ~  k^ exp 


(11) 
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where  k  is  the  integer  non-dimensional  wavenumber,  and  kg  is  the  most  energetic  wavenum¬ 
ber. 

As  Ristorcelli  and  Blaisdell^'  show,  in  a  physically  consistent  initialization  of  com¬ 
pressible  homogeneous,  isotropic  turbulence,  there  are  finite  density,  temperature,  and 
dilatational  fluctuations.  Therefore,  we  use  their  initial  conditions  and  have  compared  and 
validated  our  flow  initialization  with  that  of  Blaisdell  and  Ristorcelli.^® 

We  let  the  turbulence  evolve  to  a  realistic  state  before  allowing  the  reaction  to  progress. 
Because  the  initial  field  is  weakly  compressible,  we  use  mean  values  of  velocity  derivative 
skewness.  Si,  to  predict  the  onset  of  realistic  turbulence.  Experimental  data  presented  by 
Tavoularis  et  al}^  and  Orszag  et  al?^  show  that  Si  is  in  the  range  [—0.6  :  —0.35]  for  Rex 
near  50.  With  this  criterion,  the  turbulence  has  developed  to  a  realistic,  isotropic  state 
at  tjrt  =  0.35.  The  reaction  is  started  at  this  time,  and  we  measure  the  evolution  of  the 
turbulence  relative  to  this  time. 


We  run  cases  with  combinations  of  iSh°  =  —0.5,  1.0,  1.5,  and  2.0;  Da  =  0.5,  1.0,  and 
2.0;  Mt  =  0.173,  0.346,  and  0.519;  and  Rex  =  25.0,  34.5,  and  50.0.  We  use  nitrogen  atoms 
for  both  species  Si  and  S2,  with  initial  mass  fractions  of  csi  =  1  and  cs2  =  0.  We  choose 
Keg  =  1,  so  that  the  chemical  equilibrium  is  reached  when  csi  =  cs2  =  0.5.  kg  is  set  to 
4.  To  simplify  the  comparison  of  different  simulations,  we  freeze  jj,  at  its  initial  value,  and 
we  use  Pr  =  0.72.  We  specify  the  initial  density  as  pg  =  1.0 kg/m®.  To  be  consistent  with 
the  physical  N2  dissociation  rate,  we  use  9  =  113, 200  K. 


Having  introduced  the  modified  reaction  rate  (9),  we  must  specify  a  initial  temperature 
that  is  representative  of  a  hypersonic  boundary  layer  and  consistent  with  the  physically 
correct  reaction  rate  (5).  Using  the  law  of  mass  action  given  by  (4)  and  assuming  that  we 
initialize  the  flow  with  pure  Si,  we  can  obtain  an  expression  for  the  initial  source  term  and 
combine  it  with  (8)  to  derive  the  initial  rate  of  reaction 


l^M^-DaRTZf 

3  f^ref  R^X 


where  we  have  taken  p  from  a  power  law  (p  —  prtf{T /Tref)^)-,  Pref  from  a  Sutherland 
law,^^  and  R  is  the  universal  gas  constant.  In  order  to  obtain  realistic  conditions,  we 
choose  the  parameters  in  the  modified  reaction  rate  expression,  (9),  to  be  representative  of 
the  nitrogen  dissociation  reaction.  To  do  so,  we  use  the  N2  dissociation  rate  expression^^ 


kf  =  3.71  X  10^ 


,-113,200/T 


’/mole 


W'e  then  set  this  expression  equal  to  the  reaction  rate  in  terms  of  the  non-dimensional 
parameters.  (12),  and  find  the  temperature  that  satisfies  this  relation.  For  Mt  =  0.3, 


Da  =  0.5,  Tref  =  300K,  Rex  =  50,  and  firef  =  2.106  x  10“"^  kg/ms,  and  nitrogen 
(n  =  0.7722),  we  obtain  an  initial  value  of  T  =  9470  K.  This  initial  temperature  is 
characteristic  of  a  high-temperature  hypersonic  boundary  layer. 

VI.  Results 

We  first  focus  on  how  the  chemical  reactions  affect  the  gas  dynamics,  leaving  the 
discussion  of  chemistry  and  turbulence  interaction  for  the  next  section.  During  a  typical 
reacting  simulation,  the  average  concentration  of  species  Si  decreases  as  the  average  con¬ 
centration  of  species  S2  increases.  It  should  be  noted  that  as  the  chemical  species  reach 
the  equilibrium  mass  fraction,  csi  =  cs2  =  0.5,  the  rate  at  which  energy  is  added  to  or 
removed  from  the  fluid  approaches  zero. 


A.  Relative  Heat  Release 

In  this  section,  we  illustrate  the  effect  of  variations  of  the  relative  heat  release  at 
moderate  Da,  Mt,  and  Rex  {Da  =  1,  Mt  =  0.346,  and  Rex  =  34.5).  In  all  cases  the  non¬ 
reacting  simulation  is  shown  for  comparison  purposes.  As  mentioned  earlier,  the  relative 
heat  release  is  proportional  to  the  ratio  of  energy  released  or  absorbed  in  the  formation  of 
product  species.  Therefore,  an  increase  in  Ah°  increases  the  energy  in  the  flow  field.  This 
is  illustrated  in  Fig.  2a,  which  plots  the  temporal  evolution  of  the  average  temperature. 
Figure  2b  shows  the  temporal  evolution  of  the  turbulent  kinetic  energy  for  the  different 
values  of  Ah°.  When  Ah°  =  1,  there  is  an  initial  period  of  time  when  the  turbulence 
is  enhanced.  A  further  increase  of  the  heat  release  considerably  enhances  the  turbulent 
kinetic  energy,  maintaining  and  feeding  the  turbulence  for  a  longer  period  of  time.  For 
Ah°  —  2,  the  turbulent  kinetic  energy  reaches  high  values.  However,  at  about  t/vt  =  1.3, 
it  decays  slightly  more  rapidly  than  in  the  other  three  simulations. 

Turbulent  motion  is  dissipative,  and  needs  a  supply  of  energy  to  keep  from  decay¬ 
ing.  Moyal^^  decomposed  the  velocity  field  in  Fourier  space,  u,  into  the  incompressible 
(divergence-free),  u  ,  and  compressible  (curl-free),  u  ,  components 


A 

u  - 

It  =  u  —  u  . 


(14) 


The  dissipation  for  isotropic  turbulence  can  then  be  written  as 


e  =  8x1/ 


(F^(A;)-b  |F^(A:))F  dit. 


8 


where  and  represent  the  incompressible  and  compressible  components  of  the  shell- 
averaged  energy  spectrum 


E^  = 


E^  = 


u  •  u 
2  k  ' 

u  •  u 
2  k 


(16) 


From  (15)  we  see  that  as  the  compressibility  of  the  flow  increases,  the  turbulence  becomes 
more  dissipative.  This  explains  the  increased  rate  of  turbulent  kinetic  energy  decay  for  the 
higher  heat  release  flows. 

Figure  2c  shows  the  evolution  of  the  relative  compressible  energy  ratio 


jE^{k)dk 
J  {EC (k)  +  E^{k))  dk 


Because  the  simulations  begin  with  weakly  compressible  conditions,  when  the  turbulence 
reaches  a  realistic  state,  the  initial  values  of  x  very  low.  For  the  endothermic  reaction,  x 
is  very  small  throughout  the  simulation,  and  as  is  increased,  the  relative  compressible 
kinetic  energy  increases.  For  a  relative  heat  release  of  Ah°  =  2,  %  increases  rapidly  and 
reaches  about  20%  of  the  total  kinetic  energy  at  =  1.3.  Also  note  that  there  are 
time-dependent  oscillations  in  x  that  increase  with  the  heat  release. 

Consider  the  energy  spectrum  of  the  flow  field  decomposed  into  its  incompressible  and 
compressible  components  at  several  different  times  during  the  simulations.  Figure  3a  plots 
these  spectra  for  the  non-reacting  simulation.  We  observe  that  the  compressible  modes  are 
about  two  orders  of  magnitude  less  energetic  than  the  incompressible  modes  at  all  but  the 
smallest  scales.  However,  note  that  for  kr]  >  1^  the  length  scales  are  not  resolved  by  the 
simulation;  thus  the  increase  in  compressible  energy  spectrum  at  large  wavenumber  (small 
scales)  is  due  to  aliasing  errors.  (However,  the  level  of  error  is  small.)  As  time  evolves, 
the  compressible  energy  spectrum  decays  slightly  at  all  scales,  whereas  the  incompressible 
modes  decrease  at  the  large  scales,  and  increase  at  the  small  scales. 

Figure  3b  plots  the  energy  spectra  for  the  endothermic  simulation.  The  compressible 
modes  are  still  tw^o  orders  of  magnitude  below  the  incompressible  modes.  As  time  evolves, 
they  decrease  further  over  all  except  the  largest  scales.  Figure  3c  plots  the  energy  spectra 
for  =2.  At  all  scales,  the  energy  in  the  compressible  modes  is  about  two  orders  of 
magnitude  higher  in  this  case  than  in  the  non-reacting  case,  even  surpassing  the  incom¬ 
pressible  modes  at  the  largest  and  smallest  scales.  We  see  that  the  compressible  modes 
become  more  energetic  as  A/z®  is  increased,  but  the  incompressible  modes  are  essentially 
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unaffected  by  the  level  of  heat  release.  This  is  shown  in  Fig.  4,  which  plots  the  energy 
spectra  at  t/rj  =  1.5  for  the  non- reacting  simulation  and  the  case  with  the  highest  heat 
release.  The  incompressible  energy  spectra  are  nearly  identical.  It  is  also  interesting  to 
note  that  the  compressible  energy  spectrum  is  increased  by  the  same  factor  at  all  scales 
(except  at  small  scales  where  the  aliasing  error  is  important).  This  is  consistent  with  the 
argument  that  because  reactions  are  inherently  a  molecular-scale  process,  they  must  be 
scale  independent.^^ 

In  Fig.  5,  we  observe  a  large  increase  in  the  RMS  magnitude  of  the  temperature  fluc¬ 
tuations  when  the  heat  release  is  increased.  Recall  that  from  (9),  a  positive  temperature 
fluctuation  causes  an  exponential  increase  in  the  reaction  rate.  If  the  reaction  is  exothermic, 
heat  is  released  and  the  local  temperature  increases.  This  further  increases  the  tempera¬ 
ture  fluctuations,  and  the  process  feeds  upon  itself.  However,  because  there  is  turbulent 
motion,  the  heated  fluid  may  move  to  a  different  location  before  the  reaction  progresses 
further,  reducing  or  eliminating  the  feedback  process.  Thus,  the  interaction  between  the 
chemical  heat  release  and  the  turbulent  motion  should  depend  on  the  amount  of  heat  re¬ 
leased,  A/i°,  and  the  rate  at  which  it  is  released.  Da,.  Note  that  as  the  reaction  progresses 
toward  completion,  the  rate  of  formation  of  the  products  decreases.  Therefore,  the  net  rate 
of  heat  release  decreases  as  well,  and  there  is  less  energy  to  drive  the  fluctuations,  causing 
them  to  decrease  in  time.  The  opposite  occurs  for  the  endothermic  case:  T^j^g  decreases 
initially,  and  then  increases  in  time  as  the  turbulent  kinetic  energy  decays. 

B.  Damkohler  Number 

The  Damkohler  number  is  the  ratio  of  the  chemical  reaction  rate  to  the  rate  of  tur¬ 
bulent  motion.  For  larger  values  of  Da,  chemical  equilibrium  is  reached  in  a  shorter 
time,  which  causes  more  rapid  heat  addition  or  removal.  To  illustrate  the  effect  of  the 
Damkohler  number,  we  consider  a  case  of  moderate  heat  release,  turbulent  Mach  number, 
and  Reynolds  number  (A/i°  =  1,  Mt  =  0.346,  and  Re\  =  34.5).  Note  that  the  total 
amount  of  energy  added  to  the  flow  is  the  same  for  all  simulations. 

Consider  Fig.  6a,  which  shows  the  temporal  evolution  of  the  average  concentration  of 
the  reactant  species.  We  observe  that  as  Da  is  increased,  the  reactants  approach  equilib¬ 
rium  at  earlier  times.  Figure  6b  plots  the  evolution  of  the  turbulent  kinetic  energy.  Early 
in  the  simulations,  the  turbulent  motion  is  enhanced  with  increased  Da.  Accordingly,  the 
relative  compressible  kinetic  energy  of  the  flow  increases  with  increased  Da.  As  before, 
the  compressible  energy  spectrum  also  becomes  more  energetic,  whereas  the  incompress¬ 
ible  spectra  are  nearly  identical  for  all  simulations.  Figure  6c  illustrates  that  increased 
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Da  also  increases  the  RMS  temperature  fluctuations.  This  occurs  because  larger  reaction 
rate  cause  the  heat  released  to  remain  confined  to  a  smaller  volume  of  fluid,  enhancing  the 
subsequent  effect  on  the  turbulence.  Thus,  the  localization  of  the  heat  release  is  important 
in  maintaining  the  feedback  process. 

C.  Turbulent  Mach  Number 

To  describe  the  effect  of  the  turbulent  Mach  number,  we  choose  a  case  of  moderate 
Damkohler  number,  relative  heat  release,  and  Reynolds  number  {Da  =  1,  A/i°  =  1,  and 
Re\  =  34.5).  For  all  cases,  increasing  the  initial  Mach  number  causes  higher  values  of  (T) 
since  the  higher  turbulent  kinetic  energy  is  dissipated  and  converted  into  internal  energy. 
In  the  non-reacting  simulations,  the  relative  turbulent  kinetic  energy,  decay  rate 

depends  weakly  on  the  initial  Mf.  For  the  exothermic  simulations,  the  turbulent  kinetic 
energy  increases  initially.  This  increase  is  proportional  to  the  initial  Mf. 

Figure  7  plots  the  temporal  evolution  of  RMS  temperature  fluctuations,  T^^g,  for  the 
case  with  A/i°  =  1.  This  plot  shows  that  the  level  of  the  temperature  fluctuations  is 
proportional  to  this  is  also  the  case  in  the  non-reacting  flows.  However,  relative  to 
the  non-reacting  flow,  the  heat  release  increases  by  about  a  factor  of  two  for  all 

Mf.  Thus,  the  heat  release  amplifies  the  kinetic  energy  fluctuations,  resulting  in  higher 
temperature  fluctuations. 

D.  Reynolds  Number 

From  the  direct  numerical  simulations  we  conclude  that  there  is  very  little  effect  on  the 
decay  of  reacting  turbulence  due  to  initial  Reynolds  numbers  over  the  narrow  range  tested 
(25  <  Rt\  <  50).  We  observe  no  variation  in  the  rate  of  reaction.  At  higher  Reynolds 
numbers  the  dissipation  rate  is  reduced,  and  therefore,  the  temperature  increases  more 
slowly  and  the  RMS  temperature  fluctuations  are  somewhat  larger.  However,  these  effects 
are  insignificant  compared  to  the  variation  of  the  other  parameters. 

VII.  Interaction  and  feedback  mechanism 

The  effects  of  the  parameters  that  govern  the  interaction  between  the  turbulent  mo¬ 
tion  and  the  chemical  reactions  have  been  described  in  the  previous  section.  We  find 
that:  The  magnitude  of  the  temperature  fluctuations  increases  dramatically  with  either 
increased  exothermicity  or  increased  rate  of  reaction.  The  incompressible  energy  modes 
are  virtually  unaffected  by  the  chemical  reactions,  whereas  the  compressible  modes  are 
strongly  enhanced.  The  levels  of  compressibility  increase  with  increased  exothermicity 
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and  with  increased  reaction  rate.  From  these  results  we  can  conclude  that  there  is  a 
feedback  mechanism  between  the  turbulent  motion  and  the  chemical  reactions. 

The  velocity  decomposition  (14)  is  unique,  and  the  vorticity  is  only  a  function  of  the 
incompressible  velocity  component.  Thus,  if  the  chemical  reactions  affect  the  incompress¬ 
ible  component  of  velocity,  the  vorticity  would  also  be  affected.  Therefore,  let  us  consider 
the  evolution  of  the  vorticity  during  the  simulations.  We  can  write  the  total  change  of 
vorticity  as 

—  =  (u;  •  V)u  - w(V  •  u) -1- VT  X  V5,  (18) 

where  S  is  the  entropy  and  we  have  neglected  the  diffusive  terms.  These  terms  represent  the 
production  of  vorticity  due  to  vortex  line  stretching,  compressibility,  and  thermodynamic 
changes,  respectively.  We  can  decompose  the  third  term  by  noting  that  for  a  mixture  of 
thermally  perfect  gases 

5  =  c,lnT-i;inp+5]c,^S^,  (19) 

where  Cp  is  the  specific  heat  at  constant  pressure  of  the  gas  mixture,  R  is  the  gas  constant, 
and  ij.°  is  the  chemical  potential  of  species  s.  The  first  two  terms  form  part  of  the  entropy 
even  when  the  gas  is  not  chemically  reacting,  and  the  third  term  is  the  entropy  due  to 
chemical  reactions.  Thus,  we  can  separate  the  production  of  vorticity  due  to  thermody¬ 
namic  changes  into  two  components,  and  refer  to  them  as  the  non-reacting  entropy  and  the 
reacting  entropy.  Note  that  had  we  used  a  more  complicated  reaction,  it  would  be  more 
difficult  to  make  this  decomposition  because  the  mixture  Cp  and  R  w*ould  also  change. 

Figure  8a  plots  the  budget  of  vorticity  for  a  non-reacting  simulation  and  an  exothermic 
simulation  with  Ah°  =  2.  There  is  a  large  difference  in  the  order  of  magnitude  between  the 
terms.  The  main  contribution  to  vorticity  production  is  due  to  the  vortex  stretching  and 
turning  term  for  both  simulations.  Whereas  the  compressibility  term  remains  nearly  un¬ 
changed  for  the  non-reacting  case,  it  reaches  the  same  order  of  magnitude  as  the  stretching 
term  for  the  exothermic  simulation.  The  thermodynamic  effects  due  to  the  non-reacting 
entropy  are  larger  for  the  exothermic  simulation,  although  they  are  three  orders  of  mag¬ 
nitude  below  the  terms  mentioned  above.  The  thermodynamic  contribution  due  to  the 
reacting  entropy  is  initially  zero  for  the  exothermic  simulation,  but  increases  to  a  value 
that  is  about  an  order  of  magnitude  lower  than  the  stretching  and  compressible  terms. 
Thus,  the  main  contributions  to  the  vorticity  budget  are  the  vortex  line  stretching  and 
turning  and  compressibility  effects. 

Figure  8b  illustrates  the  same  vorticity  budget  for  the  non-reacting  simulation  and  an 
endothermic  case  with  /\h°  =  —  1.  Again,  the  main  contribution  to  vorticity  production 
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is  from  the  vortex  line  stretching  and  turning  term.  For  the  endothermic  simulation,  the 
compressible  destruction  increases  and  the  reacting  entropy  production  is  non-zero,  but 
small. 

Figure  9  plots  the  temporal  evolution  of  the  average  vorticity  magnitude  for  the  non¬ 
reacting  case  and  the  two  previous  reacting  cases.  We  see  that  the  vorticity  increases 
with  increasing  heat  release.  This  is  consistent  with  the  increased  production  mechanisms 
shown  in  Fig.  8. 

Our  explanation  of  these  results  is  as  follows.  The  heat  released  does  not  produce 
vorticity  by  itself.  Rather,  the  chemical  reactions  affect  the  turbulence  locally,  enhancing 
compressions  and  expansions.  These  compression  and  expansion  regions  interact  with 
one  another  and  generate  vorticity.  Greenberg^"*'  deduced  that  due  to  the  symmetry  of 
the  problem,  the  third  term  in  (18)  does  not  contribute  appreciably  to  the  production  of 
vorticity.  Therefore,  the  indirect  contribution  to  vorticity  from  the  chemical  heat  release 
must  be  in  the  vortex  stretching  and  turning  terms  and  in  the  production  of  vorticity  due 
to  compressibility.  Our  results  agree  with  Greenberg’s  analysis. 

The  results  also  corroborate  the  theoretical  work  of  Eschenroeder,^^  who  found  that 
the  turbulent  motion  is  fed  from  the  external  energy  source  provided  by  the  chemical 
reactions.  Our  results  show  that  the  chemical  energy  produces  local  high  temperature  and 
pressure  regions,  which  then  produce  increeised  kinetic  energy  through  the  reversible  w'ork 
term,  pV  •  u,  in  the  kinetic  energy  equation.  We  can  see  the  effect  of  this  process  on  the 
Reynolds  stress  by  considering  the  Reynolds  stress  budget 

Rij,t  =  —{Qi,j  +  Qj,i)  +  ^Rij,kk  —  Tijk,k  —  k'eij  -f  Hij.  (20) 

Where  Rij  =  ujuj  is  the  Reynolds  stress  tensor,  Qi  =  pui/p  is  the  pressure- velocity 
correlation,  Tijk  =  Uiujuk  is  the  triple  velocity  correlation  tensor,  eij  —  2ui^k'Uj,k  is  the 
homogeneous  dissipation  rate  tensor,  and  Ilij  =  p{uij  -t-  Uj,i)/ p  is  the  pressure-strain  term. 
The  pressure  transport  and  the  turbulent  transport  are  statistically  zero.  Figure  10  shows 
the  budget  of  the  R22,t  component  for  the  non-reacting  and  exothermic  simulations.  The 
pressure-strain  term  is  nearly  zero  in  the  non-reacting  case,  but  it  is  the  dominant  term  in 
the  exothermic  simulation.  Note  that  in  this  case,  there  are  oscillations  in  the  pressure- 
strain  production  term,  whose  normalized  period  is  approximately  equal  to  the  local  Mach 
number.  This  indicates  that  these  oscillations  are  a  result  of  acoustic  (pressure)  waves. 
This  further  illustrates  that  the  Reynolds  stress  generation  is  mainly  caused  by  compress¬ 
ibility  effects  in  the  reacting  case. 
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VIII.  Conclusion 

In  this  paper  we  investigated  the  interaction  of  isotropic  turbulence  and  a  chemical 
reaction  at  conditions  typical  of  a  hypersonic  boundary  layer.  We  used  a  modified  reaction 
rate  that  preserves  the  important  features  of  the  Arrhenius  rate  expression,  and  makes 
it  possible  to  perform  a  systematic  analysis.  The  flow  is  parameterized  by  four  non- 
dimensional  numbers:  the  turbulent  Mach  number,  Reynolds  number,  relative  heat  release, 
and  relative  reaction  rate. 

The  direct  numerical  simulations  show  that  the  heat  released  by  an  exothermic  reac¬ 
tion  produces  regions  of  high  temperature  and  pressure.  This  causes  localized  expansions 
which  increase  the  compressible  turbulent  kinetic  energy.  This  results  in  increased  rates 
of  kinetic  energy  decay,  vorticity  production,  and  Reynolds  stress  production.  Simulta¬ 
neously,  in  regions  of  high  temperature  the  reaction  rate  increases,  further  increasing  the 
heat  release  and  the  temperature.  Thus,  there  is  a  positive  feedback  between  the  chemi¬ 
cal  reaction  and  the  turbulent  motion.  The  strength  of  the  feedback  process  depends  on 
the  amount  of  heat  released,  the  reaction  rate,  and  more  weakly  on  the  turbulent  Mach 
number.  When  the  reaction  is  endothermic,  the  feedback  is  negative  and  temperature 
fluctuations  are  damped. 
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Figure  Captions 


Figure  1.  Time  evolution  of  the  average  N2  mass  fraction  for  simulations  of  recombin¬ 
ing  nitrogen  atoms  at  initial  Mt  =  0.519,  Rex  =  34.5,  and  Ah°  =  33,  16,  13,  and  9, 
corresponding  to  initial  To  =  2000,  4000,  5000,  and  7000  K  respectively. 

Figure  2.  Time  evolution  of  (a)  average  temperature;  (b)  turbulent  kinetic  energy;  and 
(c)  relative  compressible  kinetic  energy  showing  the  eifect  of  A.h°  for  Da  =  1,  Mt  =  0.346, 
and  Rex  —  34.5. 

Figure  3.  Three-dimensional  incompressible,  ,  and  compressible,  ,  shell- averaged 
energy  spectra  for  (a)  non-reacting;  (b)  endothermic  (A/i°  =  ~1);  and  (c)  exothermic 
(Ah°  =  2)  simulations  for  Da  =  1,  Mt  =  0.346,  and  Rex  =  34.5. 

Figure  4.  Three-dimensional  incompressible,  E^ ,  and  compressible,  E^ ,  shell-averaged 
energy  spectra  for  non- reacting  and  exothermic  (A/i°  =  2)  simulations  for  Da  =  1,  Mt  = 
0..346,  and  Rex  =  34.5  at  tjTt  =  1-5. 


Figure  5.  Time  evolution  of  RMS  temperature  fluctuations  showing  the  effect  of  A/i°  for 
Da  =  1,  Mt  =  0.346,  and  Rex  —  34.5. 


Figure  6.  Time  evolution  of  (a)  average  mass  fraction  of  reactants;  (b)  turbulent  kinetic 
energy;  and  (c)  RMS  temperature  fluctuations  showing  the  effect  of  Da  for  Ah°  =  1, 
Mt  =  0..346,  and  Rex  =  34.5. 


Figure  7.  Time  evolution  of  the  RMS  temperature  fluctuations  showing  the  effect  of  Mt 
for  Ah°  =  1,  Da  —  1,  and  Rex  =  34.5. 

Figure  8.  Normalized  vorticity,  L0tTt/u!o,  budget  for  non-reacting  (empty  symbols)  and 
(a)  exothermic  (A/i°  =  2)  (fllled  symbols);  and  (b)  endothermic  (Ah°  =  1)  (filled  sym¬ 

bols)  simulations  for  Da  =  1,  Mt  =  0.346,  and  Rex  =  34.5  ;  (□),  stretching  production, 
(A),  compressible  destruction;  (O))  non-reacting  entropy  production;  (o),  reacting  entropy 
production. 


Figure  9.  Time  evolution  of  vorticity  magnitude  showing  the  effect  of  Ah°  for  M*  —  0.346, 
Rex  =  34.5,  and  Da  =  1. 
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Figure  10.  Normalized  Reynolds  stress,  i?22,i/eo,  budget  for  non-reacting  (empty  sym¬ 
bols)  and  exothermic  (filled  symbols)  simulations  for  Mj  =  0.346,  Re\  =  .34.5,  Da  =  1, 
and  Ah°  =  2;  (o),  pressure-strain;  (v),  viscous  diffusion;  (A),  dissipation. 
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A  direct  numerical  simulation  database  is  used  to  develop  a  model  of  subgrid-scale  temper¬ 
ature  fluctuations  for  use  in  large-eddy  simulations  of  turbulent,  reacting  hypersonic  flows. 
In  these  flows,  the  reaction  rate  varies  exponentially  with  temperature,  and  thus  small 
temperature  fluctuations  may  cause  large  fluctuations  in  the  rate  of  product  formation. 
The  proposed  model  uses  a  probability  density  representation  of  the  temperature  fluctua¬ 
tions.  The  DNS  database  reveals  a  physically  consistent  relation  between  the  resolved-scale 
flow  conditions  that  may  be  used  to  predict  the  standard  deviation  of  the  Gaussian  pdf. 
The  model  is  calibrated  and  tested  by  comparison  to  simulations  of  decaying  isotropic 
turbulence.  The  conditional  single-variable  pdf  model  is  found  to  capture  the  fluctuations 
in  temperature  and  product  formation. 

Introduction 

The  boundary  layer  on  proposed  air-breathing  hypersonic  cruise  vehicles  wiU  be  tur¬ 
bulent  and  chemically  reacting.  To  aid  the  design  of  such  vehicles,  an  accurate  cfd  tool  is 
required.  Currently,  the  computer  cost  dissuades  the  usage  of  direct  numerical  simulations 
(dns)  of  hypersonic  flows  at  realistic  Reynolds  numbers,  at  the  same  time  that  it  encour¬ 
ages  the  use  of  large-eddy  simulations  (les).  Therefore,  subgrid-scales  (sGs)  models  need 
to  be  devised  to  close  the  turbulence-chemistry  interaction  in  les  of  hypersonic  flows. 

Much  work  has  been  done  in  modeling  turbulent  reacting  flows. However,  most 
of  this  work  has  been  focused  on  combustion  related  applications,  where  fluctuations  in 
the  chemical  source  term  are  mainly  caused  by  fluctuations  in  the  species  mass  fraction. 
In  hypersonic  flows,  the  equilibrium  composition  of  reacting  air  depends  strongly  on  the 
temperature.  Furthermore,  the  reactions  are  temperature  limited  and  the  reaction  rate 
depends  exponentially  on  temperature.  With  the  very  high  energies  present  in  hypersonic 
flows,  the  temperature  fluctuations  will  be  very  large,  and  result  in  large  increases  in 


1 


the  reaction  rates.  Therefore,  in  hypersonic  flows  the  temperature  fluctuations  must  be 
modeled  to  obtain  accurate  reaction  rates  and  product  formation. 

In  LES,  it  is  generally  assumed  that  the  small  scales  are  more  isotropic  and  more  uni¬ 
versal  in  character  for  different  flows  than  the  large  scales.  Therefore,  it  is  customary  to 
develop  SGS  models  in  isotropic  turbulence,  for  later  use  in  the  more  complex  flows.  In 
our  previous  work,^^  we  studied  the  interaction  between  decaying  isotropic  turbulence  and 
flnite-rate  chemical  reactions  at  conditions  typical  of  a  hypersonic  boundary  layer.  The 
results  confirmed  that  the  interaction  is  characterized  by  the  increased  or  decreased  magni¬ 
tude  of  the  temperature  fluctuations  for  exothermic  or  endothermic  reactions,  respectively. 

In  this  paper,  we  use  the  dns  database  of  isotropic  reacting  turbulence  to  develop 
and  evaluate  a  model  for  the  temperature  fluctuations.  Our  approach  builds  on  the  re¬ 
action  rate  modeling  for  combustion  applications.  In  particular,  Gaffney  et  al.^  assume 
a  probability  density  function  (pdf)  for  the  temperature  fluctuations  and  investigate  how 
the  fluctuations  affect  the  combustion  process.  The  proposed  model  has  a  similar  form, 
however  we  use  the  results  of  dns  to  calibrate  the  model.  Therefore,  the  parameters  that 
describe  the  pdf  are  based  on  detailed  turbulence  data.  This  model  represents  the  unre¬ 
solved  (subgrid-scale)  fluctuations  for  use  in  large-eddy  simulations  of  turbulence,  where 
the  small  scales  are  assumed  to  be  isotropic.  As  reactions,  we  use  either  dissociating  nitro¬ 
gen  molecules  or  dissociating  oxygen  molecules.  These  are  primary  reactions  that  occur  in 
hypersonic  boundary  layers. 

In  the  remainder  of  the  paper,  we  review  the  governing  equations,  we  Favre  average 
these  equations  and  show  that  the  unresolved  temperature  fluctuations  must  be  modeled 
for  large-eddy  simulations.  We  then  discuss  the  data  analysis  for  the  prediction  of  the 
temperature  fluctuations,  and  evaluate  the  results  by  comparison  with  the  dns  results. 

Governing  equations 

The  equations  describing  the  unsteady  motion  of  a  reacting  flow  with  no  contribution 
of  vibrational  modes  are  given  by  the  species  mass,  mass-averaged  momentum,  and  total 
energy  conservation  equations 

dps  ,  d  f  ,  _ 

dpui  d  f  ,  \ 


(1) 


where  ps  is  the  density  of  species  s;  Ws  represents  the  rate  of  production  of  species  s  due 
to  chemical  reactions;  Ps  is  the  density  of  species  5;  uj  is  the  mass-averaged  velocity  in  the 
j  direction;  Vsj  is  the  diffusion  velocity  of  species  s;  p  is  the  sum  over  species  s  density;  p 
is  the  pressure;  nj  is  the  shear  stress  tensor  given  by  a  linear  stress-strain  relationship;  Qj 
is  the  heat  flux  due  to  temperature  gradients;  hs  is  the  speciflc  enthalpy  of  species  s;  and 
E  is  the  total  energy  per  unit  volume  given  by 

E  =  '^  PsCvsT  +  ^puiUi  +  '^psh°,  (2) 

S  S 

where  c„s  is  the  specific  heat  at  constant  volume  of  species  s;  and  /i°  represents  the  heat 
of  formation  of  species  s. 

We  consider  a  simple  reaction  system:  dissociating  nitrogen  molecules.  This  is  one 
of  the  primary  reactions  that  occur  in  hypersonic  boundary  layers,  but  more  importantly, 
it  serves  as  a  good  example  for  the  form  of  the  chemical  source  term.  The  dissociation 
reaction  may  be  written  symbolically  as: 

N2  +  M  ^  2N  -}-  M,  (3) 


where  M  represents  a  collision  partner,  which  is  either  N2  or  N  in  this  case.  The  source 
terms  for  N2  and  N  can  be  written  using  the  law  of  mass  action 
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and  wn,  =  -wn.  kf  and  ki,  are  the  forward  and  backward  reaction  rates  respectively; 
These  are  written  as 

kf=CT^e-^^^, 

k  - 

where  C,  77,  and  6  are  experimentally-determined  constants, and  K^q  is  the  known 
temperature-dependent  equilibrium  constant. 


For  a  two-species  mixture,  the  diffusion  velocity  can  be  accurately  represented  using 
Pick’s  law 

(6) 


D^Cs 


where  Cs  =  Ps/p  is  the  mass  fraction,  and  D  is  the  diffusion  coefficient  given  in  terms  of 
the  Lewis  number 


Le  = 


pDPr 

P 


(7) 
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where  Pr  is  the  Prandtl  number,  ^  is  the  viscosity,  and  Le  is  taken  to  be  unity,  so  that 
the  energy  transport  due  to  mass  diffusion  is  equal  to  the  energy  transport  due  to  thermal 
conduction. 


Favre  average  of  conservation  equations 

To  discuss  how  the  conservation  equations  must  be  modeled,  we  need  to  separate  the 
resolvable  fluctuations  from  those  that  the  grid  cannot  resolve.  In  compressible  flows,  it 
is  convinient  to  use  Favre-filtering  to  avoid  the  introduction  of  subgrid-scale  terms  in  the 
equation  of  conservation  of  mass.  A  Favre-filtered  variable  is  defined  as: 


f  =  ^ 

P  ' 


(8) 


where  the  over-bar  denotes  a  filtering  operation  that  maintains  only  the  large-scales.  Ac¬ 
cordingly,  the  Favre-filtered  equations  describing  the  unsteady  motion  of  a  reacting  flow 
with  no  contribution  of  vibrational  modes  are  given  by  the  species  mass,  momentum,  and 
total  energy  equations: 

—  (  nr  .\  -I-  - 

dt 


dt 


^  [pUiUj  +  p5ij  -I-  Tij  -  aij^  =  0  (9) 

d 

+  ^ ((-^  +  +  pj  +  K^sj^j  -  0 


where  Cs  =  Ps/p  is  the  Favre-averaged  mass  fraction,  and  the  SGS  diffusion  velocity  has 
been  introduced 

/  \ 

(10) 


Vsj  —  pi^s^j 

In  the  momentum  equation,  rij  is  the  SGS  stress 

X-ij  —  p^^LiXtj 

and  has  been  simplified  according  to 


<Tj  =  P: 


dui 


+  P' 


_  duj  2  _  duk 


-  7;P: 


Sij  , 


(11) 


(12) 


dxj  '  dxi  3  dxk 

thus,  the  viscosity  is  assumed  to  be  a  weak  function  of  the  mixture  mole  fractions,  and  its 
fluctuations  are  neglected.  In  the  energy  equation,  we  have  introduced  the  average  ratio 
of  specific  heats,  7,  the  SGS  heat  flux 


Qj  ~  p{^h'^  > 


(13) 
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and  the  SGS  kinetic  energy  flux 


Sj  —  2  Uk^k^j^  • 


(14) 


These  ..terms  are  derived  by  considering  the  Favre  average  of  the  energy- velocity  product 


ExLj  —  QC-usCsXijT  p'^s^ 

=  Euj  -1-  T,j  +  T,sh°Vsj  -1-  p'EsCvs  (csTuj  -  Cgfuj'j 
Euj  +  ^Qj  +  Sj  +  Y^sh^Vsj 


(15) 


where  we  have  assumed  that  the  mass  fraction  fluctuations  are  not  directly  correlated  with 
the  temperature-velocity  variations.  The  energy  gives  the  temperature 


E  —  p^T  -(-  pUiUi  -f"  pTg/igCg  -[-  pE^Cy^^g  ) 

z  ^ 


(16) 


where  =  Tcg  -  fcg.  And  the  pressure  is  given  by 


p  =  pRT  -I-  pRUs^s 


(17) 


The  results  from  DNS  of  isotropic  turbulence  show  that  the  terms  involving  the  SGS 
temperature-species  correlation,  $s,  are  negligible  in  comparison  to  the  rest  of  the  terms 
in  the  equations  of  state  and  total  energy.  Thus,  if  this  result  prevails  in  the  case  of  the 
hypersonic  boundary  layer,  will  drop  out  in  the  large  eddy  simulations. 

Finally  we  must  consider  the  chemical  source  term,  Ws,  in  the  continuity  equation. 
The  strong  temperature-dependence  of  the  source  term  can  be  seen  by  representing  the 
variables  as  their  Favre-averaged  and  subgrid-scale  values:  T  =  T and  Cg  =  Cs-f  c®^®. 
The  source  term  (4)  to  first  order  in  the  subgrid-scale  quantities  is: 


Wfi^o  =  W/N2  +  '*^/N2  (^(;h  +  ^) 


WbN2  =  WbN2  +  'f^6N2 


T 

T  ^  T 


j^sgs 

~T~ 

J'sgs 


(18) 


where  we  have  made  use  of  the  fact  that  and  we  have  assumed  that  the 

variation  of  in  temperature  is  relatively  weak,  which  is  a  good  approximation  at  very 
high  temperatures.  The  variations  of  caused  by  temperature  fluctuations  may  be 
especially  large  because  9  is  typically  an  order  of  magnitude  larger  than  the  temperature. 
Thus,  to  obtain  accurate  Wg  we  must  model  the  subgrid-scale  temperature  fluctuations. 
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This  derivation  shows  that  the  following  terms  must  be  modeled: 


Vsj  =  SGS  diffusion  of  species  s, 

Tij  =  303  stress, 

Qj  =  SGS  heat  flux, 

T,j  =  SGS  kinetic  energy  flux, 
j-sfifs  _  gQg  temperature  fluctuations. 


(19) 


Much  work  has  been  done  in  modeling  nj  and  whereas  very  little  effort  has  been 

given  to  modeling  Vgj,  Ej,  and  However,  Gao  and  O’Brien^^  briefly  discuss  the 

modeling  of  Vsj  in  the  context  of  a  dynamic  SGS  model,  and  Reveillon  and  Vervisch^® 
discuss  the  modeling  of  Ej .  The  current  paper  includes  the  development  of  a  model  for 


Simulations  of  reacting  turbulence 


In  this  section  we  present  the  numerical  simulations  for  three-dimensional,  compress¬ 
ible,  homogeneous,  isotropic,  reacting  turbulent  flow.  We  use  the  same  numerical  method 
and  isotropic  initialization  as  Martin  and  Candler.  In  particular,  we  use  a  sixth-order 
accurate  finite  difference  scheme, and  a  fourth-order  accurate  Runge-Kutta  time  integra¬ 
tion.  The  simulations  are  performed  on  grids  with  96^.  The  fluctuating  fields  are  initialized 
using  Ristorcelli  and  BlaisdelR^  consistent  initial  conditions. 

The  non-dimensional  parameters'^  governing  the  turbulence-chemistry  interaction  are 
the  turbulent  Reynolds  number,  the  Mach  number,  the  Damkohler  number,  the  relative 
heat  release  and  the  equilibrium  constant,  namely 


Rex 


pu'X 

IJ- 


Ah°  =  - 


CyT  + 


(20) 


and  is  a  non-dimensional  function  of  the  temperature  and  determines  the  equilibrium 
mixture  composition,  q  is  the  RMS  magnitude  of  the  fluctuation  velocity;  a  is  the  speed 
of  sound;  n'  is  the  RMS  turbulent  velocity  fluctuation  in  one  direction;  A  is  the  Taylor 
microscale;  r*  is  the  turbulent  time  scale;  Tc  is  the  chemical  time  scale;  and  Ah°  is  the 
heat  of  the  reaction. 


We  choose  initial  Re\  =  34.5,  (p)  =  0.5  kg/m^,  and  (T)  =  3000  K  for  N2  dissociation, 
and  (T)  =  2000  K  for  O2  dissociation.  At  these  initial  temperatures  the  reactions  are 


6 


exothermic.  We  run  simulations  with  combinations  of  initial  Mt  —  0.35,  and  0.52,  and  N2 
or  O2  mass  fractions  in  [0:1].  The  various  initial  mass  fraction  conditions  lead  to  distinct 
evolutions  of  the  nonequilibrium  chemistry.  Therefore,  different  turbulent  flow  fields  are 
obtained.  The  values  of  M°,  Da,  K^q,  and  Mt  vary  physically  during  the  simulations. 


Prediction  of  temperature  fluctuations 


In  this  section,  we  describe  the  development  of  a  model  for  the  subgrid-scale  temper¬ 
ature  fluctuations.  We  represent  the  temperature  fluctuations  as  a  pdf.^^  During  a  les  we 
would  obtain  the  temperature  fluctuations  by  sampling  from  the  pdf.  From  dns,  we  use 
histograms  and  statistical  quantities  to  obtain  the  shape  of  the  pdf  for  the  temperature 
fluctuations.  This  determines  the  number  of  parameters  needed  to  predict  the  distribution. 
We  analyze  the  data  from  the  dns  by  plotting  the  parameters  that  describe  the  pdf  against 
Mt,  Da,  Ah°,  and  Keq  evaluated  using  ensemble- averaged  quantities  during  the  simula¬ 
tions.  Using  the  bivariate  scatterplots,  we  find  a  functional  relation  for  the  parameters 
describing  the  distribution. 

Figure  1  plots  the  skewness,  St',  and  flatness,  Ft',  factors  for  typical  reacting  sim¬ 
ulations.  After  the  initial  transient  St'  is  nearly  zero,  and  Ft'  is  about  3;  these  values 
indicate  that  the  temperature  fluctuations  may  be  characterized  by  a  Gaussian  distribu¬ 
tion.  Two  parameters  describe  a  Gaussian,  the  mean  and  the  standard  deviation.  We 
choose  to  look  for  a  distribution  of  the  normalized  temperature  fluctuations,  T' /{T).  In 
this  way,  the  mean  is  zero  and  we  only  need  to  predict  the  standard  deviation,  which  is 
given  by  T^m^/{T).  When  analyzing  the  results,  we  only  consider  the  data  after  the  initial 
transient,  which  corresponds  to  t/xt  >  0.25. 

The  scatterplots  indicate  that  T^ms/{T)  changes  considerably  with  Da  and  Mt  during 
the  simulations.  is  a  function  of  temperature  and  also  affects  the  standard  deviation. 
Figure  2  shows  a  functional  relation  of  T^ms/(^)  ^^e  nitrogen  dissociation  simulations. 

We  observe  that  when  plotting  T’^ms/(^)  versus  the  relative  heat  release  times  a  length 
ratio,  X/Ie,  the  data  collapses  into  a  power  fit.  The  length  ratio  is  defined  from  the 
following  relation  between  the  governing  parameters: 


X/Ie^ 


u'xt 

ar* 


—  Mt  Da\  log(/C 


eq) 


(21) 


u'  Tt  is  the  distance.  A,  traveled  by  a  particle  of  fluid  moving  at  the  speed  of  the  turbulent 
intensity.  W^e  have  redefined  the  chemical  time  scale  (20)  to  take  into  account  the  effect 
of  the  equilibrium  constant,  r*  =  rd  log(A'eg)l-^  Then,  a  r*  is  the  distance,  /e,  covered 
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by  the  acoustic  radiation  due  to  temperature  fluctuations,  and  gives  the  deflnition  of  the 
expansion  length,  Ie-  Moin^'^  argues  that  the  chemical  process  must  be  independent  of  the 
turbulence  length  scale,  since  it  occurs  at  a  molecular- level.  Martin  and  Candler^ ^  find 
that  the  turbulence-chemistry  interaction  is  independent  of  the  Reynolds  number  for  the 
range  of  Rex  that  they  consider.  Thus,  it  may  be  surprising  that  the  length  ratio  gives 
a  relation  for  T’^ms/(^)-  However,  the  Taylor  microscale  is  defined  by  the  velocity  scale 
of  the  inertial  range  eddies  and  the  time  scale  associated  with  the  dissipative  eddies.^® 
Therefore,  the  Taylor  microscale  is  not  associated  with  a  particular  eddy  size.  The  dns 
results  show  that  \/Ie  is  always  less  than  1.  As  X/Ie  approaches  one,  T^^^/{T)  is  large 
which  indicates  a  strong  turbulence-chemistry  interaction.  This  situation  corresponds  to 
the  case  where  the  fluid  motion  covers  a  similar  distance  to  that  covered  by  the  acoustic 
radiation  induced  by  the  temperature  fluctuations.  However,  as  X/Ie  approaches  zero  the 
thermal  expansion  outruns  the  fluid  motion  and  the  interaction  is  weak.  This  is  reflected 
in  the  small  values  of  T^^s/ (^)-  addition,  T^^s/ (^)  affected  by  the  amount  of  relative 
heat  released  to  the  flow.^^  Thus,  the  length  ratio  is  modulated  by  Ah°  to  give  the  relation 
for  the  standard  deviation: 

Kujm  =  A(M^XIlEf,  (22) 

where  A  and  B  are  constants  listed  in  Table  I. 

While  the  above  relation  (22)  and  physical  argument  hold  for  exothermic  reactions, 
the  treatment  of  endothermic  reactions  must  be  a  different  one.  When  the  flow  reaches  the 
required  internal  energy  to  permit  the  dissociation  of  molecules  the  reactions  become  en¬ 
dothermic  and  the  flow  becomes  incompressible.’^^  Thus,  the  magnitude  of  the  temperature 
fluctuations  is  reduced  by  the  chemical  reactions  and  the  acoustic  radiation  is  nullified. 
In  this  situation,  the  temperature  fluctuations  are  proportional  to  M^.  This  is  also  re¬ 
flected  in  Fig.  2,  where  the  filled  symbols  follow  a  straight  vertical  line  and  correspond 
to  the  part  of  the  simulation  where  the  reactions  are  endothermic.  Consider  Fig.  3  which 
plots  T^ms/(^)  versus  for  nitrogen  reactions  at  different  initial  Mach  numbers.  The 
straight  line  indicates  the  values  of  the  standard  deviation  when  the  chemical  reactions 
are  endothermic.  And  the  empty  symbols  that  lie  outside  the  straight  line  correspond  to 
^RMs/ (^)  when  the  chemical  reactions  are  exothermic.  Thus  for  endothermic  reactions  the 
standard  deviation  is  given  by: 

TLms/{T)  =  constant  Mf,  (23) 

where  the  constant  of  proportionality  is  0.9  for  nitrogen  dissociation,  and  0.48  for  oxygen 
dissociation.  Figures  4a  and  4b  show  the  curve  fits  for  the  oxygen  dissociation. 
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A  simple  single- variable  pdf  representation  of  T'  gives  no  distinction  between  positive 
or  negative  temperature  fluctuations.  Thus,  given  the  pdf  for  T' / (T),  we  use  a  conditional 
sampling  to  obtain  the  temperature  fluctuations  at  each  point  of  the  flow.  Figure  5  shows 
scatter  plots  of  the  temperature  fluctuations  against  the  pressure  fluctuations,  P' / (P)  at 
t/Tt  =  0.5  and  2.  We  observe  that  T'  is  positively  correlated  with  P'  for  most  part  of  the 
flow  domain  during  the  dns.  In  the  case  of  a  LES,  we  can  extrapolate  the  sign  of  T'  using 
a  double  filter^'^  to  obtain  the  sign  of  the  resolved  P' . 


Model  evaluation 

We  consider  the  situation  where  the  resolved-scale  temperature  is  known  from  the 
solution  of  the  conservation  equations,  but  the  temperature  fluctuations  are  unknown.  At 
every  time  step  during  an  exothermic  simulation,  we  compute  using  the  corre¬ 

sponding  power  law  in  terms  of  Mt,  Da,  Ah°,  and  Kgq  evaluated  using  ensemble-averaged 
quantities.  At  every  point  in  space,  we  conditionally  sample  from  the  Gaussian  pdf  de¬ 
scribed  by  r^Ms/(^)  mean.  Based  on  the  new  T' ,  we  compute  the  new  tempera¬ 

ture.  This  procedure  results  in  a  physically  correct  statistical  average  of  the  subgrid-scale 
temperature  fluctuations  and  the  chemical  source  term.  In  the  following  paragraphs  we 
evaluate  the  model  by  comparison  to  the  dns  results. 

For  nitrogen  dissociation  with  initial  Mt  =  0.52  and  initial  N2  mass  fraction  of  0.3  and 
0.8,  Fig.  6a  plots  the  time  evolution  of  T^^s-  Let  us  first  consider  the  dns  results.  For  both 
simulations  the  reactions  are  initially  exothermic  and  the  temperature  increases.  During 
the  simulation  with  initial  Na  mass  fraction  of  0.3,  the  temperature  reaches  an  optimal 
value  for  the  recombination  of  nitrogen  atoms,  and  the  reaction  becomes  endothermic  near 
f/rt  =  1.  Thus,  reaches  nearly  325  K  in  less  than  0.25  t/rj.  But  decays  rapidly 
during  the  rest  of  the  simulation,  showing  weak  turbulence-chemistry  interaction  at  later 
times.  In  contrast,  for  the  simulation  with  initial  mass  fraction  of  0.8,  the  reaction  is 
still  exothermic  passed  t/rt  =  2.  reaches  roughly  175  K  in  0.25  and  is  maintained 

during  the  simulation.  Consider  now  the  predictions  given  by  the  temperature  fluctuation 
model.  For  both  simulations  the  model  predicts  T'^us  accurately.  For  (cn2)o  =  0.3  the 
expression  for  the  standard  deviation  is  switched  from  (22)  to  (23)  near  f/rt  =  0.5.  When 
Ah°X/lE  <  Mf  Figure  6b  shows  that  the  model  represents  acceptably  during  the 
simulation.  It  should  be  noted  that  the  magnitudes  of  given  by  the  dns  are  very 
small.  Such  small  changes  in  species  fluctuations  do  not  contribute  significantly  to  the 
production  of  species.  And  therefore  are  not  important.  Figure  6c  shows  that  the  model 
represents  correctly  during  the  simuliition. 
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To  further  illustrate  the  performance  of  the  model,  we  consider  two  cases  for  the 
dissociation  of  O2.  Figure  7  plots  the  modeled  and  dns  results  for  simulations  with  initial 
Mf  =  0.52  and  initial  O2  mass  fraction  of  0.3  and  0.8. 

For  the  previous  simulations,  Fig.  8  plots  the  correction  factor  appearing  in  (18) 
and  given  by  the  dns.  This  Fig.  shows  that  the  temperature  fluctuations  contribute 
significantly  to  the  species  production  early  dnring  the  simulations,  when  the  reactions 
are  exothermic.  However  when  the  reactions  become  endothermic,  such  are  the  cases  with 
(cn2)o  =  (coa)©  =  0.3,  the  temperature  fluctuations  are  small  and  do  not  contribute  to  the 
production  of  species  considerably. 

The  model  also  shows  good  agreement  with  other  turbulent  quantities  such  us  ensem¬ 
ble  and  RMS  values  of  T,  Ws,  Cs,  and  turbulent  kinetic  energy  decay.  Higher  moments  such 
as  dilatation  values  are  not  well  predicted.  When  the  model  is  turned  on,  the  dilatation 
values  over  shoot.  This  is  a  direct  conseciuence  of  the  procedure  that  we  use  to  test  the 
model. 


Conclusions 

A  PDF  approach  is  used  to  develop  a  subgrid-scale  temperature  fluctuation  model  for 
the  closure  of  the  chemical  source  term  in  large-eddy  simulations  of  hypersonic  flows.  We 
use  a  DNS  database  to  obtain  the  shape  of  the  distribution  for  T'/(T).  We  find  that  the 
distribution  is  Gaussian  and  can  be  described  by  the  standard  deviation,  T^^Ms/iT),  and 
zero  mean.  For  exothermic  reactions  the  magnitude  of  the  temperature  fluctuations  is 
enhanced  and  the  standard  deviation  can  be  expressed  as  A{/S.h°\/l£)^ .  Where  \/Ie  is 
directly  obtained  from  the  non-dimensional  governing  parameters,  and  it  represents  the 
ratio  of  the  local  characteristic  distance  covered  by  the  fluid  motion  to  the  characteristic 
distance  covered  by  the  acoustic  radiation.  For  endothermic  reactions,  the  flow  becomes 
more  incompressible  and  the  temperature  fluctuations  scale  with  .  We  find  that  T' 
is  positively  correlated  with  P'  for  most  part  of  the  flow  domain.  We  use  this  condition 
to  obtain  the  sign  of  the  modeled  temperature  fluctuations.  We  evaluate  the  model  by 
comparison  to  dns  and  find  good  agreement  of  results.  Previous  results^^  indicate  that 
the  chemistry-turbulence  interaction  is  independent  of  the  turbulent  length  scales.  Thus, 
we  expect  the  constants  appearing  in  (22)  and  (23)  to  be  independent  of  the  grid  size  filter. 

In  the  context  of  a  hypersonic  boundary  layer,  the  values  of  the  governing  parameters 
and  therefore  the  length  ratio  can  be  obtained  by  averaging  along  path  lines  or  planes 
parallel  to  the  wall.  Then,  the  temperature  fluctuations  may  be  obtained  by  conditionally 
sampling  from  the  probability  density  function  described  by  the  calibrated  power  curves 
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for  exotliorinic  reactions.  For  endothermic  reactions,  the  magnitude  of  the  temperature 
fluctuations  is  small,  and  will  not  contribute  significantly  to  the  production  of  species. 
Therefore,  the  temperature  fluctuations  may  not  need  modehng  in  a  situation  where  the 
reactions  are  endothermic.  It  still  remains  to  be  determined  if  a  Gaussian  probability 
density  function  can  be  used  to  accurately  represent  the  temperature  fluctuations  in  a 
boundary  layer.  Owen  et  al.^®  find  that  for  a  cold  wall  boundary  layer  at  Moo  —  tb® 
probability  density  function  for  the  temperature  fluctuations  is  skewed  near  the  wall. 
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Table  Captions 


Table  I.  List  of  constants  for  T^ms/(^)  =  A  (A/i°  X/Ie)^. 
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TABLE  I 


Na 

Oa 

Mt 

A 

B 

A 

B 

0.35 

0.016 

0.199 

0.062 

0.179 

0.52 

0.043 

0.225 

0.048 

0.306 
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Figure  Captions 


Figure  1.  Time  evolution  of  temperature  fluctuations  skewness,  St',  and  flatness,  Ft', 
factors  for  N2  dissociation  at  initial  Mt  =  0.52. 

Figure  2.  Standard  deviation  versus  the  modulated  length  ratio,  Ah°\/lE,  for  N2  disso¬ 
ciation  at  initial  Mt  =  0.35, 0.52. 

Figure  3.  Standard  deviation  versus  for  N2  dissociation  at  initial  Mt  =  0.17,0.35 
and,  0.52. 

Figure  4.  Standard  deviation  for  O2  dissociation  (a)  versus  the  modulated  length  ratio, 
Ah°X/lE,  at  initial  Mt  =  0.35, 0.52;  and  (b)  versus  M^  at  initial  Mt  —  0.17, 0.35  and,  0.52. 

Figure  5.  Scatter  plot  of  the  normalized  temperature  fluctuations,  T' /{T),  and  the 
normalized  pressure  fluctuations,  P'/{P),  for  N2  dissociation  at  initial  Mt  —  0.52,  and 
(cn2)o  =  0.8,  after  (a)  t/rt  =  0.5;  and  (b)  t/rt  =  2. 

Figure  6.  Model  performance  for  N2  dissociation  at  initial  Mt  =  0.52  and  (cn2)o  =0.3, 
0.8.  Time  evolutions  of  (a)  rms  temperature  fluctuation;  (b)  rms  of  species  fluctuations; 
and  (c)  ensemble  average  of  the  source  term,  {w]^^)/{wi<[^)o. 

Figure  7.  Model  performance  for  O2  dissociation  at  initial  Mt  =  0.52  and  (coj)©  =0.3, 
0.8.  Time  evolutions  of  (a)  rms  temperature  fluctuation;  and  (b)  ensemble  average  of  the 
source  term,  {wo^)/{wo^)o. 

Figure  8.  Time  evolution  of  factor  (-^  appearing  in  (18)  for  N2  dissociation 

with  (a)  (cn^)©  =  0.3,  solid;  (b)  (CN2)©  =  0.8,  dashed;  and  for  O2  dissociation  with  (c) 
(002)0  =  0.3,  long-dashed;  and  (d)  (coa)©  =  0.8,  dashed-double  dotted;  all  simulations 
start  with  initial  Mt  =  0.52. 
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